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Abstract 

A generalization of the Renormalization Group, which describes order-parameter 
fluctuations in finite systems, is developed in the specific context of percolation. 
This "Stochastic Renormalization Group" (SRG) expresses statistical self-similarity 
through a non-stationary branching process. The SRG provides a theoretical basis 
for analytical or numerical approximations, both at and away from criticality, when- 
ever the correlation length is much larger than the lattice spacing (regardless of 
the system size). For example, the SRG predicts order-parameter distributions and 
finite-size scaling functions for the complete crossover between phases. For percola- 
tion, the simplest SRG describes structural quantities conditional on spanning, such 
as the total cluster mass or the minimum chemical distance between two bound- 
aries. In these cases, the Central Limit Theorem (for independent random variables) 
holds at the stable, off-critical fixed points, while a "Fractal Central Limit Theo- 
rem" (describing long-range correlations) holds at the unstable, critical fixed point. 
This first part of a series of articles explains these basic concepts and a general 
theory of crossover. Subsequent parts will focus on limit theorems and comparisons 
of small-cell SRG approximations with simulation results. 
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1 Introduction 



The theory of critical phenomena is mostly concerned with phase transitions in 
the thermodynamic limit of infinite systems. In this limit, an order parameter 
has unambiguous values in each phase and is singular at a precise critical 
point of the control variable (temperature, concentration, etc.) [1-4]. In finite 
systems, however, there are two basic modifications to the thermodynamic 
picture: (i) The critical point becomes blurred since all statistical averages 
are analytic functions of the control variable, and (ii) the order parameter 
fluctuates around its mean value. Accordingly, the order parameter in a finite 
system is described by a probability distribution, which depends non-trivially 
on the control variable and the system size. 

As arguably the simplest critical phenomenon, percolation [5-7] provides a 
natural setting within which study order-parameter fluctuations in finite sys- 
tems. Moreover, in recent years there has been growing practical interest in the 
distribution functions of structural quantities, such as the minimum or maxi- 
mum chemical distance between two terminals [8-17], the masses of spanning 
clusters [11-13,18], and the mass of the largest cluster [19-22]. These ran- 
dom quantities not only act as order parameters for the phase transition, but 
also play a crucial role in both dynamical [23,24,16] and mechanical [19,21] 
properties of heterogeneous materials [25] . 

There is no existing theory to predict order-parameter statistics for arbitrary 
values of concentration and system size, but considerable progress has been 
made in understanding the limiting behavior in each thermodynamic phase. 
Classical limit theorems for independent random variables hold in the normal 
phases, where correlations are negligible (by definition). In the supercritical 
phase, the masses of large fragments of the infinite cluster [26-28] and other 
structural quantities [29,30] satisfy the Central Limit Theorem, i.e. the dis- 
tribution is Gaussian, and the mean and variance both grow linearly with 
the system size. In the subcritical phase, the mass of the largest cluster is 
governed by the theory of extremes of independent random variables [19,20]. 
In this case, the distribution has the Fisher- Tippett (or Gumbel) form with 
non-universal corrections due to discrete-lattice effects, while the mean grows 
logarithmically with system size (with a bounded variance) [20] . 

Computer simulations of percolation have shown that various order parame- 
ters in the critical phase have different universal distributions with stretched- 
exponential tails (when scaled to have unit mean) [8,9,11,18,14-16,22]. These 
signatures of what we call the "Fractal Central Limit Theorem" [31] can be 
predicted by analyzing stationary branching processes [32-34], as was recently 
pointed out by Hovi and Aharony [13]. Aside from an early analysis of mean- 
field percolation on the Gayley tree [10], these authors made the first (and. 
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to the author's knowledge, only) systematic use of branching processes to ap- 
proximate critical limiting distributions, although they focused mainly on the 
tails of spanning-cluster mass distributions [13]. 

In this first part of a series of articles, we present the "Stochastic Renormaliza- 
tion Group" (SRG), a theoretical formalism based on non-stationary branch- 
ing processes which describes the crossover of order-parameter distributions 
between different thermodynamic phases. The SRG expresses the statistical 
self-similarity of finite systems in which the correlation length, ^, and linear 
system size, L, are both much larger than the lattice spacing, a. The critical 
^ L) and off-critical ^ L) fixed points of the SRG correspond to various 
limit theorems for the order parameter. In the latter case, small correlation 
lengths ~ a) are also allowed. Although the basic idea of the SRG applies 
to critical phenomena and random fractals in general, we develop it here on 
the specific context of percolation. 

The article is organized as follows. In section 2, we set the stage by review- 
ing classical renormalization-group (RG) concepts in percolation. In section 3, 
we present the SRG for any statistical quantity conditional on spanning the 
system, such as the total mass of spanning clusters or the minimum chemical 
distance. In section 4, we show how the SRG generally describes crossover. 
In sections 5 and 7, we briefiy explain why these the off-critical and critical 
fixed points are associated with the CLT and the FCLT, respectively. In sec- 
tion 8, we summarize the SRG "fiow" in the space probability distributions 
for percolation quantities dependent upon spanning, and we conclude by dis- 
cussing the general relevance of the SRG for fractals and phase transitions 
in finite systems. Subsequent parts will analyze the fixed points in detail, in- 
cluding more comparisons with simulation data. The subtle case of the largest 
cluster (which may not span), the most natural order parameter in finite sys- 
tems [19,20,35-37], will also be addressed. 



2 RG Concepts in Percolation 

2.1 Self- Similarity at the Critical Point 

The basic idea behind all RG methods for critical phenomena is the statistical 
self-similarity of large systems near a critical point [2-4] . This essentially fol- 
lows from dimensional analysis [38]: In the limit a/L ^ and a/^ 0, there 
are only two (out of originally three) relevant length scales, L and ^, so all 
systems with the same dimensionless ratio, L/^, must have approximately the 
same statistical properties (upon a suitable normalization). The self-similarity 
becomes exact in the thermodynamic limit, L — > oo, as the critical point is 
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approached, ^ — > oo, with the scahng variable, L/^, held constant. 

The notion of self-similarity is exploited by the RG to calculate nontrivial 
properties of the critical point, which cannot be obtained by simple dimen- 
sional analysis. For example, in the case of percolation, according to the ar- 
guments above, a large system of N — mn sites must have approximately 
the same connectivity statistics as a coarse-grained system of n ^ 1 large 
"cells" (or "blocks") containing m 3> 1 sites each, provided that an effective 
occupation probability for cells, = Rm{p), is defined so as to preserve the 
dimensionless ratio, L/^: 

where d is the system dimension {N — L"^). This is analogous to Kadanoff's 
original block renormalization of the coupling constant for magnetic spins [4] . 

If the exact dependence of the correlation length on p were known, then the 
cell occupation probability would be given by 

«»w = r'(||i). (2) 

The trick is to now reverse the logic: Prom a suitable definition of the cell 
occupation probability, an estimate of the correlation length, can be 

obtained by successively "renormalizing" the occupation probability, 

Pmn — RniPm)-! (3) 

starting from pi — p. By construction, an appropriate choice of Rnip) yields a 
valid asymptotic approximation in the large-cell limit, \\m.n^oo i^^\p) = i{p)i 
although it is not trivial to make this choice for a given problem such that the 
convergence is fast enough to allow the use of small cells. 



2.2 Asymptotic Approximations 

Estimates of correlation length and p^ itself are obtained from the RG, Eq. (3), 
as follows. As shown in Fig. 1, for any suitable choice of Rn{p), the recursion 
(3) converges to one of three possible fixed points in the thermodynamic limit, 
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Fig. 1. Typical trajectories of the renormalized cell occupation probability pN given 
by Eq. (3) for the four-site {N = 4^) RG scheme described below in Eq. (11) starting 

at the critical fixed point, pi = pi^^ = (\/5 — l)/2 (solid line), a nearby subcritical 
value, p'l = p[^^ — 10^^ (dotted line), and a nearby supercritical value, pi = p^^ + 10~^ 
(dashed line). An associated estimate of the correlation length, ^^^^(p'l), in the sub- 
critical case is labeled as the system size iV = ^'^ where the renormalized occupation 
probability pN drops below 0.5. This is very close to a correlation length estimate, 
^^^\pi), for the supercritical case where pat rises above 0.7. 

7 oo: 



lim = < 



in) 
c 



for < p < p 

pin) fQj, p _ pin) 

1 for < p < 1 



(4) 



the stable subcritical fixed point, p = 0; the stable supercritical fixed point, 
p = 1; or the unstable critical fixed point, p — p^^\ which is an 7i-cell estimate 
of the true infinite-system value of Pc- As long as the RG transformation, is 
chosen to satisfy 



lim Rnip) 



for < p < Pc 
c for p = Pc 

1 for Pc < p < 1 



(5) 



for some constant < c < 1, then the RG approximation converges in the 
large-cell limit, limn^ooPi"-* = Pc- As shown in Fig. 1, for p ^ Pc''\ an estimate 
of the correlation length, ^'■"•'(p), can be defined as the (interpolated) linear 
system size, n-'/'^, where the cell occupation probability, p„j , first passes within 
a given tolerance of a stable fixed point as j — > oo. 
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The RG introduces a useful artificial length scale, the linear cell size, b {n = b'^). 
In order for the lattice-spacing a to not affect the RG calculation of critical 
behavior, we must have b ^ a, but the underlying assumption of self- similarity 
also requires b <^ C{p), which is guaranteed near the critical point, where 
CiPc) — CO- General considerations of dimensional analysis imply that such 
a divergence must have the form of a power law, ^{p) ^ {p — Pc)'" (e.g. 
for p p^), but dimensional analysis alone cannot determine the critical 
exponent u. The RG, however, easily produces an asymptotic estimate, 
for a given choice oin — b'''. 

Substituting an approximation of the correlation-length singularity, 

(6) 



into Eq. (1) with m — n yields 



Rn{p)-P^:'' -n'^^'^p-pf^), 



(7) 



from which we identify, n^'^'^'"' = The resulting RG estimate of the 

critical exponent is 



(Note that, since Rn{p) is generally a polynomial, it is analytic at p^"\ which 
implies that the correlation length diverges with the same exponent for p — > 
p^ and p ^ Pc ■) Such RG estimates of critical exponents are believed to 
converge to their exact values in the large-cell limit, provided that Eq. (5) 
holds, consistent with our analysis below and extensive numerical evidence in 
the literature. 



2.3 The Renormalization Group 



Algebraically, the recursion in Eq. (3) represents a group of operators acting on 
the occupation probability (more precisely, a semigroup without an inverse.) 
There are at least two equivalent ways to construct this "renormalization 
group". First, by using the same number of cells at each level of recursion, the 
renormahzed occupation probabilities, p„j = TZ^f^p, for the sequence of system 
sizes, N = n, n^, n^, n'^, . . ., can be generated by iterates of the function, Rn{p), 
acting on the site occupation probability, p e [0, 1]. In this case, the group 



6 



operators, 



nf^ ^R°J = RnOR^O... 



oR, 



(9) 



(j compositions of Rnip)) are mappings from [0, 1] to [0, 1] which depend on 
the cell size, n. 

Alternatively, one could view Rnip) as the initial condition for a group of 
operators acting on the space of functions, {/ : [0, 1] [0, 1]}, by setting the 
number of cells to equal the number of sites at each iteration. For example, 
the renormalized probabilities, p^2j = T^jRnip), for the sequence of system 
sizes, N — n,n'^,n'^,n^, . . ., can be generated by group operators defined as 
follows. 



Although the initial condition depends on the function Rn{p), in this case the 
RG is always the same. As such, broad sets of initial conditions can yield the 
same limiting behavior, and in this way the RG explains, at least qualitatively, 
the universality of critical exponents. 

2.4 Cell Occupation by Spanning 

These abstract concepts would not be terribly important were it not for the 
fact that they can be put to practical use by clever definitions of the cell 
occupation probability. In the spirit of this commemorative issue, we focus 
on the seminal contributions of Reynolds, Stanley, and Klein [39-41]. Before 
their work, a few simple, small-cell decimation procedures were proposed [42- 
45], but these authors first defined a cell to be "occupied" it is spanned by a 
cluster from one side to another in at least one direction. This choice captures 
the essence of percolation as a phase transition in connectivity. Moreover, it 
should always produce convergent RG approximations because Eq. (5) holds 
rigorously for any spanning rule that, as n — > 00, requires a cluster to extend 
across an infinite distance [6]. 

For a given choice of spanning rule, the spanning probability, Rnip), can be 
evaluated exactly for small cells or approximately by Monte-Carlo simulation 
for large cells [39-41]. For example, for a square site or bond lattice in two 
dimensions, the function Rn{p) could be the probability that a cluster connects 
opposite boundaries of a square cell in one direction, both directions, or either 
direction. Once Rn{p) is obtained, estimates of Pc and u are given by the 
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expressions above. The fractal dimension Dfoi incipient infinite clusters can 
also be estimated by generalizing the RG to inchide "ghost bonds" between 
non-nearest neighbor sites, although below we derive a different estimate which 
does not require any ghost bonds. 

For illustration purposes, throughout this paper we will consider the following, 
very simple example. For site percolation on the square lattice, a reasonable 
RG for a small cell of only n = 2 x 2 sites is based on the probability of 
spanning in one direction [40,41], 

Ra{p) = + Vg + 2p^g^ (11) 



(where q — 1 —p) which is the sum of probabilities to have spanning clusters of 

size four, three, or two, respectively. Solving the quartic equation, Ri{p) — p, 
yields (curiously) the golden mean as an estimate of the critical point, 

p(4) = ^ ^ 0.618034, (12) 



which is fairly close to the best available numerical value [46], Pc — 0.59274621(13). 
The critical exponent estimate from Eq. (8), 



.(4) 



1.63529, 



(13) 



is not as close to the presumably exact value, u = 4/3, although larger cells 
slowly improve this estimate [41]. Of course, more sophisticated and accurate 
RG schemes exist for various percolation problems, but this simple example 
illustrates the power of RG methods to approximate nontrivial quantities, such 
as the critical exponents and p,. itself, with remarkable ease. Using the same 
example, we will produce RG estimates below for many other quantities, both 
at and away from the critical point, with little additional effort. 

In the limit of large cells, n — > oo, numerical evidence suggests that RG esti- 
mates converge to their exact values (whenever known), and the accuracy of 
small-cell approximations is often quite remarkable [39,41,47,13,59]. In spite of 
these successes, however, the validity of the large-cell RG based on spanning 
has recently been debated in the context of critical spanning probabilities [48- 
56,47,57-59]. Strictly speaking. Equation (3) seems to predict that the critical 
spanning probability in the infinite-system limit is precisely equal to pc- On 
the other hand, simulations on the square site lattice support another theo- 
retical prediction (from conformal field theory [60]) that the critical spanning 
probability in two dimensions is exactly 1/2, which is somewhat less than p^. 
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There is no contradiction, however, since the RG only predicts that the span- 
ning probabihty at the n-cell estimate is equal to p^J^\_ which need not 
hold in the limit of infinite cells [47]. In fact, when the efi'ect of the finite 
cell size is properly considered as — > Pc, the RG actually predicts the 
universal spanning probability, including its subtle dependence on the system 
shape [48,60,55,57,58] and boundary conditions [53,54,58] with good preci- 
sion [47,59]. In summary, the general validity of the RG based on spanning 
is now fairly well established, and thus we use it as a starting point for our 
theory below. 



2.5 Self- Similarity Away From the Critical Point 



The RG is often invoked to motivate (but not to calculate) the "data collapse" 
of any quantity, Q{L,p) depending on system size L and concentration p onto 
finite-size scaling functions of the form, 

Q{L,p)r.L^f^(J-]^ (14) 



where the exponent, Df, (which is the fractal dimension, in the case of span- 
ning clusters) and the dimensionless scaling function, $(a;), do not depend on 
L or p. (Here and below, we always choose units to set a = 1.) As explained 
above, this result is simply a consequence of dimensional analysis. 

There is no existing theory to predict finite-size scaling functions. Although 
the fractal dimension can be calculated from the RG, the scaling function 
is usually obtained by the ad hoc fitting of numerical data [5,7,61,20]. This 
approach, however, is quite tedious and fails to take advantage of the full 
power of the RG. 

Although RG methods calculations have mostly been applied to the critical 
phase, where ^ ^ L (or ^ = oo), the general arguments given above only 
require that the correlation length be large, ^ ^ 1, but not necessarily larger 

than the linear system size L. Therefore, in principle the RG can not only 
motivate, but also predict, finite-size scaling functions, as we demonstrate by 
a suitable generalization below. 

For RG approximations using cells of size n — V^, the scaling law, Eq. (14), is 
replaced by 

«(L.p) . L-r*.") (15) 
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where ^^"^(p), Dy are estimates based on the n-cell level of coarse 

graining. The scahng function, ^^"^(x), can be determined at any value of its 
argument, L/^^"-\p), by an appropriate renormalization of p to maintain this 
ratio in the limit N = L'^ = ^ oo, as explained below. Since h/^{p) — * in 
this limit, the general arguments above imply convergence, lim^^oo ^'^"H^) — 

If the system size is large, L ^ 1, and one is not interested in accurately de- 
scribing crossover (which also requires being close to the critical point, ^ 3> 1), 
then the RG can actually produce accurate approximations far from the criti- 
cal point, even when the correlation length is small, ^(p) ~ 1. As emphasized 
in Ref. [20], if a system is partitioned into cells with b ^ then each cell 
makes independent contributions to any quantity of interest for the whole 
system, which are easily combined using standard techniques from probability 
theory. The possibility of accurate approximations in this limit is a conse- 
quence of self-similarity near the stable, off-critical fixed points of the RG, 
where ^ < 6 < L. 



3 Stochastic Renormalization Group for Spanning- Cluster Masses 

3.1 The Mathematical Statement of Statistical Self-Similarity 

A key ingredient missing in classical RG methods is the randomness of the 
order parameter in finite systems. We now reformulate the RG in probabilistic 
terms for the case of percolation, using a general formalism which could also 
be applied to other critical phenomena. Let Xjv(p) be an order parameter 
conditional on spanning, such as the total mass of spanning clusters or the 
minimum chemical distance between opposite boundaries, for a large system 
of size A?" » 1 at a given concentration p, sufficiently close to the critical point 
that discrete lattice effects are small, ^{p) ^ 1. Consider coarse-graining such 
a system of size N into n cells of size m. 

The powerful constraint of statistical self-similarity discussed above implies 
that, upon a suitable renormalization of the occupation probability, the ran- 
dom number of occupied cells, Xn{pm), has the same distribution as in a system 
of size n with the original value of p. Unlike the sites or bonds in the original 
system, however, the coarse-grained cells do not have equal masses. Instead, 
each occupied cell contributes a random mass, X^{p), sampled independently 
from the appropriate distribution for a smaller system of size m. 

This motivates a recursive construction we refer to as the "Stochastic Renor- 
malization Group" (SRG): The total spanning-cluster mass is given by a re- 



10 



cursive random sum of random variables, 



(16) 



1=1 



where X^^ denotes the ith independent sample from the mass distribution for 
each celL Random cluster masses, Xn{p), X„2(p), Xn3{p), . . ., at the n-cell 
level of approximation can be generated recursively using Eq. (16), which is 
fully determined by the cell probability distribution. 



and the cell occupation probability, Rn{p), from Eq. (3). These simple quan- 
tities can be either calculated exactly for small cells or approximated numeri- 
cally for large cells. The SRG then extrapolates the distribution to larger finite 
systems in an arbitrary stage of crossover between phases. 

The SRG approximates the statistics of cluster masses through random geo- 
metrical constructions, as illustrated in Figure 2. The coarse-graining proce- 
dure described above involves discarding contributions from unoccupied cells, 
and the figure only shows the sites and cells which survive the entire process 
(from right to left). Equivalently, clusters are created by recursively refining 
occupied cells (from left to right). In each row, the picture on the right shows 
a random cluster of sites, and the pictures to the left show the hierarchical 
structure of connectivity among the sites at different scales of coarse-graining. 
At scales smaller than the correlation length, these artificially generated clus- 
ters are random fractals characteristic of the critical phase, while at larger 
scales in normal phases they are objects of integer dimension, as explained 
below. Although the connectivity of these artificial clusters is not required 
explicitly (as shown in the figure), their mass distributions must nevertheless 
approach those of the real spanning clusters in each phase in the large-cell 
limit, purely as a consequence of statistical self-similarity. The crossover be- 
tween phases can also be described if the correlation length is much larger 
than lattice spacing. 

3.2 Analogy with Population Growth 

In mathematical terms. Equation (16) describes a stochastic branching pro- 
cess [32-34], which is notably non-stationary since the generating distribu- 
tion, fn{pm,x), is different at every stage of recursion. Since the nineteenth 
century, stationary branching processes have been analyzed as models of pop- 
ulation growth, e.g. in the contexts of human family trees [62], animal pop- 




(17) 
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Fig. 2. Three random realizations (in each row from left to right) of a nine-cell SRG 
branching process for the total cluster mass conditional on spanning in the vertical 
direction. Top row: the supercritical phase, showing the crossover to critical behavior 
at scales smaller than the correlation length. Middle row: the critical phase, showing 
complete statistical self-similarity. Bottom row: the subcritical phase, also showing 
the crossover to critical behavior. (Note that in these examples the fractal dimension 
has been reduced well below its actTial value for 2d percolation, Df = 91/48, to 
illustrate the crossover regime with greater clarity.) 

ulations [63] and nuclear chain reactions [64]. It does not appear, however, 
that any mathematical studies have addressed the particular type of non- 
stationarity in Eq. (16), which has multiple fixed points corresponding to 
different thermodynamic phases. 

The analogy with population growth provides a delightfully simple way to 
understand the SRG: The hierarchical structure of power-law correlations is 
analogous to the genealogical structure of populations. This connection is easily 
understood from the graphical constructions of Fig. 2. Each cluster of sites 
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is a graphical representation of a "family tree" . The "ancestors" of a site are 
the coarse-grained cells in which it is contained. From the pictures, it is clear 
that any random changes to the occupied cells ( "ancestors" ) at a given length 
scale (or "generation") greatly affect the number of occupied cells and sites 
("descendants") at all smaller length scales. Therefore, even though all of the 
random variables in the branching process are independent, long-range spatial 
correlations, analogous to "family ties" , are induced on the final "population" 
of sites. 

Although the analogy with population growth is conceptually very useful, the 
SRG branching process involves several new features not typically studied in 
the context of populations, and thus it requires new mathematical results. For 
example, the large-cell limit (of infinitely many "children" per individual) is of 
primary interest for the SRG. More importantly, the structure of spatial corre- 
lations would corresponds to a very unusual "family tree" where randomness 
in the number of children only appears gradually after a certain generation. 
In the SRG, this occurs when the cell size is comparable to the correlation 
length (the middle column in Fig. 2). Below this scale, spatial correlations 
induced on the sites decay as a power law, reflecting the absence of a charac- 
teristic length, while above this scale correlations decay exponentially because 
the number of "children" is no longer random. This gradual variation in the 
probability distribution for "births" in the branching process causes a phase 
transition as the linear system size passes through the correlation length, as 
explained below. 



3.3 Recursions for Distributions and Generating Functions 

Since all the random variables in the SRG branching process are independent. 
Equation (16) takes the form [32], 

fmn{P:X)= fn{PmJ)fm{P:X): (18) 

in terms of probability distribution functions, where f*^ = f * f * ■ ■ ■ * f (j 
times) denotes the jth auto-convolution, defined by f*g{x) = J2y f{x—y)g{ii). 
In general, the bounds of summation, x~ and x^, are set by the distribution 
function, /„(p, x), which also initializes the recursion in Eq. (18). In the case of 
percolation quantities subject to a spanning condition, the minimum number 
of occupied cells is of the order of the linear cell size, = 0{b), e.g. = 6 = 
^i/d hypercubic cells, and the maximum number occupied cells is always 
equal to the total number of cells, x'^ — n. 
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For both analytical and numerical calculations with the SRG, it is convenient 
to work with various generating functions [32], which are essentially discrete 
Fourier or Laplace transforms. For example, Equation (18) takes the simple 
form, 

fmn(p, Z) = fn(Pm, fm(p, z)), (19) 

in terms of probability generating functions, 

oo 

f{z)^j:f{x)z\ (20) 



x=Q 



Other useful forms of the SRG, 

/mn(s) =/m(Pm, log /m(P,s)) (21) 
fmn(p,s)^fn(Pm,fm(p,s)), (22) 

are expressed in terms of moment generating functions, 

= (23) 

or cumulant generating functions, 

/>)=log/», (24) 



respectively. These transformations reduce the complicated sum of convolu- 
tions in Eq. (18) to a simple composition of generating functions. In physics 
parlance, this amounts to changing the representation of the SRG from real 
space to various definitions of reciprocal (or momentum) space. 



3.4 The Stochastic Renormalization Group 



Let us consider the algebraic structure of the transformations involved in the 
branching process, Eq. (16). There are various ways to iterate the recursions 
above, as we now illustrate using the cumulant-generating-function represen- 
tation, Eq. (22). For example, as in Eq. (9) the same number of cells n can be 
used at each level of recursion, which yields approximations, 

fn^ (P) = fn(Pnj-i ) ° fn {Pn^-^ ) ° ■ ■ ■ ° fn{Pn) ° fn{p) , (25) 
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for the sequence of system sizes, = n, n^, n^, n^, . . .. Alternatively, as in 
Eq. (10) the number of cells can be set equal to the number of sites at each 
level of recursion, which yields the (identical) approximations 



fn^(p) ^ fn{Pn) ° fn{p) 

= fn{Pn3) O fn{Pn^) O fn{Pn) ° fn{p) 

fn^+^ (P) = (Pn^^ ) ° L^^ (P) 

= /n(]3„2J+i_l) O /n(p„2J+i_2) « • ■ • ° fn{Pn) O fn{p) 

for the sequence of system sizes, = n,n'^,n'^,n^, .... At a fixed point, the 
same set of RG operators {TZj} as in Eq. (10) now acts on initial conditions, 
fn{p, s), in the space of functions, {/ : [0, 1] x C i— >• C}. 

If p is not a fixed point, then the functional transformations are somewhat 

more complicated, but they still have the structure of an Abelian semigroup 
(hence the name "Stochastic Renormalization Group") in terms of a set of 
operators, {T^^^*^}, defined by 



• {R, /(P)} = {R o R, f(R(p)) o f(p)}. (26) 

These functional operators act on initial conditions, {Rn{p), fn{p, s)}, in the 
function space, {R : [0, 1] ^ [0, 1]} x {/ : [0, 1] x C ^ C}, 

{it:„2. (p), (p, s)} = 7^f «{i?„(p). Up, (27) 



where Rn{p) is the spanning probability and /n(p, s) is the cumulant generat- 
ing function for the spanning cluster mass in a cell of size n. 

The fact that the SRG can be expressed as a universal group acting on a 
space of initial conditions suggests that the crossover of distributions between 
different phases shares the same universality properties as the fixed points. 
Therefore, the shapes of scaling functions and the limiting distributions in the 
critical and crossover regimes should depend on embedding dimension, but not 
on the microscopic details of the lattice (although we do not claim to rigorously 
prove this universality here, even within the SRG formalism). Below, we will 
argue that boundary conditions also affect the universal limiting distributions 
and thus also the scaling functions for crossover. 
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4 A Theory of Crossover Phenomena 



4.1 Two Types of Fluctuations 



It is straightforward to derive recursions for the moments, 



I^uAp) = {XniPn = I fn{p, 1) = ( ^ J /n(P>0), (28) 



or the cumulants, 



CnAp)= Up,^\ (29) 



by differentiating the various generating-function recursions above. For exam- 
ple, the mean, //„(p) = A*n,i(p) — Cn,i(p), satisfies the recursion, 

l^mnip) = l^n{Pm)l^m{p): (30) 



which simply confirms our intuition that the mean total mass is equal to the 
mean number of occupied cells times the mean mass per cell. The variance, 
c'n(p)^ = Cn,2{p), satisfies another recursion, 

CTmnipf = IJ'n{Pm)(7m{pf + CTniPmf fJ'mipf , (31) 



which is coupled to that of the mean. 

The recursion for the variance, Eq. (31), has an important physical interpre- 
tation. The first term on the right side is the variance within each cell times 
the mean number of occupied cells. This represents the well-known additivity 
of variance for sums of a fixed number of independent random variables. The 
second term is the additional variance due to fluctuations in the number of 
occupied cells. As explained above, this kind of fluctuation induces hierarchical 
(power-law) correlations between the microscopic site variables. The gradual 
disappearance of this second source of variance at scales larger than the cor- 
relation length is precisely what determines the crossover away the critical 
phase. 

The higher cumulants are important indications of this crossover since they 
measure departures from a Gaussian distribution, so we give a few more recur- 
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sions implied by the SRG branching process. The skewness, cn,3{p), satisfies, 

(P)o-m(p)^ + Cn,3{Pm)lJ^m{pf, (32) 



and the kurtosis, cjv,4(p), satisfies, 



CmnAip) ^ l^n(Pm)CmAp) + (^n{Pmf {^l^m{p)Cm,3{p) + 3(7^ (p)^) 
+'^Cn,3{pm)lJ'm{pY(7m{pf + Cn,4(Pm)A*m(p)^- 



(33) 



On the right side of each of these recursions, the first term represents the 
additivity of the cumulants of the mass within each cell over the expected 
number of occupied cells, as in sums of independent variables. The remaining 
terms represent additional contributions due to fluctuations in the number of 
occupied cells. As wc elucidate below, the phase transition occurs as a result 
of the competition between these two kinds of fluctuations, representing short- 
range and long-range correlations, respectively. 

4-2 Crossover of Probability Distributions 

As discussed in section 2, large systems near the critical point which have the 
same ratio, L/^{p), are in the same stage of crossover, i.e. they are statistically 
identical. This situation can be conveniently arranged within the formalism of 
the SRG, where the analogous ratio is L/^^^\p): Different systems are in the 
same stage of crossover if they have the same cell occupation probability on the 
largest cells, i.e. at the highest level of coarse-graining. In particular, a sequence 
of systems of increasing size can be frozen in the same state of crossover by 
adjusting p at each stage of recursion to set the occupation probability on the 
largest cells to a given constant, < 77 < 1, which acts as the finite-size scaling 
variable. 

For example, if the cells have reached size m — for a system of size N — 
n^^^, then we adjust p to the following value 



(the inverse of Rn{p) applied j times). This choice ensures that the correlation 



p = R-^{r,) = R-'oR-^o...R-\r^) 



(34) 




L L b 
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where n — and N — L'^. 

The crossover from subcritical to critical to supercritical is controlled by vary- 
ing rj from to p^"^ to 1, respectively. For a given value of 77, if we define, 

h^+i (77, s) = {R~^ (77) , s) , (36) 
then the SRG takes the form, 

hn<V) - fniv) O fn(R-\v)) 
kMv)-fn{v) O fniK'iv)) 0...0 fn{Rn\v)) (37) 

in the cumulant-generating-function representation. By construction, the in- 
verse iteration of the occupation probability converges to the critical fixed 
point, 

limi?;^(,7)=p(") (38) 



for < 77 < 1, so the iteration for /ijv(^, s) approaches a stationary functional 
iteration, which has a well defined limit (after an appropriate rcscaling, dis- 
cussed below). If we next take the large-cell limit, n oo, then convergence 
is attained because the first term in the resulting infinite product tends to the 
true distribution inside a large cell of size 7^ — > oo, while each of the remaining 
terms tends toward a constant, /in(p^"\s), because (again, by construction) 
lim„^oo -Rn (^) — pi"^- Therefore, in the large-cell limit the SRG produces 
asymptotic approximations of scaling functions and probability distributions 
in the crossover regime. 



4-3 Finite-Size Scaling Functions 

An analytical formula for the scaling function of the kth cumulant, 

c.,M ~ (i"")'*?" (^) . (39) 
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can be derived from the n-cell SRG. The trick is to choose rj according to 
Eq. (35) and thus cast the scahng function in the equivalent form, 



1(4) ■ («' 



Although the subcritical {p < p^"^) and supercritical (p > p^"^) scahng func- 
tions are different in terms of the variable x = L/$^^"'\p), there is a single, con- 
tinuous scahng function in terms of the variable rj G [0, 1]. From the auxiliary 
function, ^^"^(?7), the usual scaling function for the A;th moment, $[,"''(2;), can 
be determined from Eq. (40), once the appropriate n-cell correlation length, 
^("■^(p), is measured (e.g. as in Fig. 1). 

The auxiliary scaling function, ^^"''(r]), is easily calculated (at least numeri- 
cally) from Eq. (37). For example, consider the case of the mean, k = 1. After 
j iterations of the SRG, an approximation of the scaling function, ^'^"■'(r/), is 
given by 



/^n(p^"V+^ /^n(pi"^) /^n(p^"^) /^n(pi"^) /^n(p^"^: 



where the rescaling factors in the denominator are chosen so that 

*(")(p(")) = $(")(0) = 1. (42) 

This is also required for convergence as j — > cx) (with 77 fixed) since R^-'iv) ~^ 
p^"\ In this limit, we obtain a formula for the n-cell SRG approximation of 
the infinite-system scaling function. 



*S"'(^) = n (43) 

=0 fJ'niPc ') 



It is straightforward to prove that this infinite product converges if //^(p^"^) > 
and i?^(p^"^) > 0, as required for the critical exponents to be finite and 
nonzero. 

By comparing Eq. (39) and Eq. (41), we have 



which yields a simple estimate of the fractal dimension, 

# = dlog„/.„(p(")). (45) 
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It is apparent that this estimate of the fractal dimension converges to the 
correct value in the large-cell limit, as long as p^^^ Pc, because it is equivalent 
to the usual finite-system estimate used in numerical simulations, iiniPc) ^ 
■nPsI^. Equation (45) has also been derived by Hovi and Aharony using similar 
methods, albeit vahd only at the critical point [13] (see also below). Here, 
however, we also predict scaling functions and probability distributions for 
the complete crossover between phases. 

As n — > oo, the first term in the infinite product of Eq. (43) converges 
to the true scaling function, and all other terms converge to unity, since 
Xmin^ao^^ij]) = Pc"'' f*^^ < ?7 < 1 according to Eq. (5). This suggests 
that the SRG approximation of the scaling function converges in the large-cell 
limit. A very similar convergence proof is carried out in detail in Part II for 
an infinite-product formula for the strength of the infinite cluster. 



4-4 A Simple Example 



To illustrate the SRG approach, we compute estimates of finite-size scaling 
functions for the total mass of clusters spanning a square system in one direc- 
tion, using the very simple four-cell scheme introduced above in section 2.4. In 
this case, the SRG is based on the conditional probability generating function, 

- , , 2p^q^z'^ + Ap^qz^ + p^z^ , 

^^^^'") = ' ^ ^ 



where R\{j)) is defined in Eq. (11), which yields the following polynomials for 
the moments, 

, , 2'=+y ^2 ^ 3^=4^^^ + 4V 
/^4,.(P) = . (47) 



Prom these expressions it is straightforward to compute (at least numerically) 
the auxihary scaling functions ^^"^(77), but it would also be necessary to cal- 
culate the correlation length, ^^^^{p) to obtain the usual scaling functions. 

In practice, it is easier to iterate the cumulant recursions (numerically) starting 
from two particular values of p very close to p^f> (one above and one below) , as 
in Pig. 1, and then rescale by appropriate values of ^,^"'\p) for the chosen values 
of p. This method effectively samples the true continuous scaling functions at 
discrete points V /i'^'^\p)-, which are very closely spaced (in log coordinates). 
Por our four-cell example, the results of this method are shown in Pig. 3. 



20 



100 



10 



|X /L 



0.1 



0.01 



0.001 



1 ' 1 
■ (a) 

- 


1 ' 1 ' 1 , 
(4) 

slope d - „ ' ' 

-- 


- 

: P>pr--- 


(4)'''-. 


- p=p:'— 


slope 1 - Df - 


P < p(4) 




1 , 1 


1 , 1 , 1 '■ 



0.1 r 



0.01 r 



o /L ' 



0.001 



0.0001 



0.00001 



0.0001 



0.0001 



0.01 



0.01 



100 



10000 



1 ' 1 ' 1 ' 1 ' 1 




^ X (4) " 




"-. slope 




(4) " V 




slope 1/2 -Dj \ 


■ P>pr--- 




; p=p':— 




p < p(4) 

1 


1 , 1 , 'i. 



100 



10000 



Fig. 3. Four-cell SRG estimates of crossover functions for (a) the mean and (b) the 
standard deviation of the total mass of clusters spanning a square system in one 
direction. 

The scaling functions in Fig. 3 exhibit a smooth crossover from one power 
law to another (straight lines on a log-log plot) as the system size exceeds 
the correlation length, which is the correct qualitative behavior. The quanti- 
tative accuracy of these particular approximations is not very good, however, 
because even the fractal dimension (limiting slope) is poorly described. For 
this example, Equation (45) predicts D^^^ ^ 1.611, for the mass of spanning 
clusters, compared to the presumably exact value, Df = 91/48 ~ 1.896. 



As with the critical exponent, u, four cells {b — 2) is not enough to reason- 
ably approximate the fractal dimension, Df, but these estimates improve very 
quickly with increasing cell size. For example. Equation (45) predicts fractal 
dimensions for various structural quantities on the square bond lattice with 
less one percent error using remarkably small cells {b < 4) [13]. Likewise, 
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preliminary numerical results with the SRG using larger cells on the square 
site lattice also suggest that the approximation of scaling functions improves 
quickly with cell size, concomitant with the improvement in describing the 
fractal dimension. This ongoing work will be reported elsewhere, but at least 
the analysis here predicts convergence in the large-cell limit (provided that 
the scahng functions actually exist). 



5 The Supercritical Fixed Point 

5.1 Asymptotic Independence 

A detailed analysis of the fixed points of the SRG will be presented in upcom- 
ing parts of this series, so here in the remaining sections we simply give the 
basic ideas. 

Above the critical point of the RG, p > p^P'K once the linear cell size exceeds 
the correlation length, m^^'^ ^ ^^""Kp), the renormalized cell occupation prob- 
ability approaches the supercritical fixed point, Pm 1- In this limit, the 
sums in Eq. (16) no longer have a random number of terms because Xn{l) is 
equal to a constant, x^, with probability one, i.e. 

fn{l,z) = z^K /n(l,s) = e^"^ and /„(l,s)=4s. (48) 

The suppression of randomness in the number of "children" of the SRG branch- 
ing process corresponds to the removal of hierarchical correlations at scales 
larger than the correlation length. In the infinite-system limit, independent 
contributions of mass from all cells are simply accumulated at each level of 
recursion. In other words, the SRG at the supercritical fixed point reduces to 
a non-random sum of independent random variables, 

i=l 

For example, at scales larger than the correlation length in the supercritical 
phase, all = n = b'^ cells make independent contributions to the total mass 
of spanning clusters, while only a straight line of = 6 cells contributes to 
the minimum chemical distance in a hypercubic system. (The former case is 
illustrated in the top row of Figure 2.) 
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Fig. 4. Supercritical and subcritical probability distributions, g{z), for the total mass 
of clusters spanning a square system in either direction (scaled to have zero mean 
and unit variance). Numerical simulation data for square- lattice site percolation is 
compared with the analytical prediction of the SRG at the off-critical fixed points, 
i.e. the Gaussian, /^/\/27r. The numerical data reflects averaging over 10 million 
samples for a system of size N = 200 x 200 for the supercritical case and 20 million 
samples of size N = 100 x 100 for the subcritical case. 



5.2 Central Limit Theorem 



It is well known that non-random sums of independent random variables obey 
the Central Limit Theorem (CLT), which is characterized by a Gaussian hm- 
iting distribution with linear scaling of the mean and variance [32]. In the 
case of the SRG, "the GLT holds" means that the distribution of the scaled 
random variable 

Xn{p) - hn{p) 
Zjv — [oo 



converges to a standard normal distribution (i.e. a Gaussian with zero mean 
and unit variance). If equality holds in Eq. (49) and thus Eq. (48) applies, then 
the CLT is straightforward to prove from the cumulant generating function 
representation, 

fmn{s,p) = X+fm{s,p), (51) 



which implies the additivity of cumulants (the crucial step in most proofs of 
the CLT). This recursion is only asymptotically valid, however, so it is not 
clear a priori that the convergence is fast enough to ensure that the CLT 
holds. 
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Nevertheless, in Part II it is shown that the CLT docs in fact hold for the super- 
critical fixed point of the SRG. This result is consistent with Newman's proofs 
of the CLT for the mass of finite fragments of the infinite cluster [26,65,28]. 
The SRG prediction is also in excellent agreement with our own new simu- 
lation results (using the methods of Ref. [20]) for the total mass of clusters 
spanning a large square system in either direction, shown in Fig. 4. 



6 The Subcritical Fixed Point 

For p < the occupation probability renormalizes to zero because the 
spanning probability tends to zero (exponentially) once the cell size exceeds 
the correlation length. In this limit, all cells become "unoccupied", and hence 
spanning clusters are extremely rare. When spanning clusters do occur, how- 
ever, they tend to include the minimum number of cells, x~, in a straight 
line across the system {x~ — b), as illustrated in Fig. 2. Indeed, this intuition 
is supported by the SRG for the conditional distribution of spanning cluster 
masses: At the subcritical fixed point, the SRG again reduces to a non-random 
sum of random variables, 

i=l 

since X„(0) = x~ with probability one. 

As in the supercritical case, this suggests that the CLT holds. The simulation 
data in Fig. 4 is also consistent with this claim, although it is difficult to gen- 
erate enough subcritical spanning clusters to attain good convergence of the 
distribution. (Much better convergence could presumably be achieved using 
a "go-with-the-winners" algorithm [66].) To the author's knowledge, the va- 
lidity of the CLT for various quantities conditional on spanning in subcritical 
percolation has not previously been reported. 

Since spanning clusters are very rare in the subcritical regime, a much better 
order parameter for the entire phase transition in finite systems is the mass 
of the largest cluster [35-37,19,20]. At the subcritical fixed point, the mass 
distribution of the largest cluster does not obey the CLT, and instead the 
Fisher- Tippett limit theorem of extreme statistics holds [19,20]. On the other 
hand, the CLT does hold for the largest-cluster mass at the supercritical fixed 
point because the largest cluster and various spanning clusters are all the same 
in the supercritical limit. To unify the present theory with these intriguing 
results revealing subcritical/supercritical asymmetry, in Part IV the SRG is 
extended to the more subtle case of the mass of the largest cluster (which may 
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not span). 



7 The Critical Fixed Point 

7.1 Complete Statistical Self- Similarity 

At the critical fixed point, p = p^\ the SRG reduces exactly to 

X^^^^X^\ (53) 

i=l 

which asserts that, at every scale, the random cluster mass in a system 
of A*" = mn sites has the same statistics as the sum of independent mass 
contributions X^ from a random number X„ out of n cells of size m. (In this 
section, we will drop the p arguments in our notation with the understanding 
that p = p'^\) As was recently first pointed out by Hovi and Aharony [13], the 
n-cell iteration of Eq. (53), with the distribution of X^ fixed in the upper limit 
of the random sums, corresponds to a stationary Watson-Galton branching 
process [67,33,34,32]. This classical discrete branching process describes the 
growth of a population of individuals, each of whom has a random number of 
children, X^, sampled independently from the same distribution. 

The self-similar iteration of Eq. (53), 

(54) 

1=1 



starting from the initial condition X„ (where A^ = n^^) corresponds to a hier- 
archical sampling of the Watson-Galton process: The distribution of grandchil- 
dren is that of the children's children, the distribution of great, great grand- 
children is that of the grandchildren's grandchildren, etc. In this form, the 
SRG at the critical fixed point clearly expresses statistical self- similarity at 
all length scales because the coarse-grained mass "on the cells" (the random 
number of terms in the sum) has the same distribution as the random mass 
contributions from each cell. 

In terms of cumulant generating functions, the SRG at the critical fixed point 
is simply a functional iteration, 

fmn — fn^ fmi (^^) 
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starting from the initial condition fn{s)- From Eq. (30), this imphes that the 
mean scales with system size, N — L'^ — , as a power law, 

1^^ = ii^^y = L^f = (n^/ 'y (56) 

where the n-cell estimate of the fractal dimension is given once again by 
Eq. (45). 

(Prom Eq. (56), we see that //„ > 1 is a necessary condition for the expected 
value of the critical order parameter to grow with system size, or equivalently, 
for the fractal dimension to be positive. In the context of population growth, 
this is the case where the expected population size grows exponentially with 
the number of generations, e.g. a supercritical nuclear reaction. Conversely, 
the case, //„ < 1, models exponentially fast extinction of a population, e.g. a 
subcritical nuclear reaction. As such, probability theorists refer to these two 
types of branching processes as "supercritical" and "subcritical" , respectively, 
and while the marginal case /^^ = 1 is naturally termed "critical" [32-34]. 
Unfortunately, this standard mathematical terminology is poorly suited here 
because the SRG at the critical fixed point is actually a "supercritical branch- 
ing process". In fact, the non-stationary SRG branching process is actually 
"supercritical" for any thermodynamic phase since the expected value of the 
order-parameter always increases with system size.) 

7.2 The Fractal Central Limit Theorem 

Although the SRG corresponds to a new kind of non-stationary branching pro- 
cess, various existing mathematical results for simple Watson-Galton branch- 
ing processes [67,33,34] have direct relevance for the critical fixed point [13]. 
For example, the "central region" of the limiting distribution is at the same 
scale as the mean. In other words, the scaled random variable, = X^/fiisf, 
has a stationary limiting distribution, 

lim Prob.(X„. <z^i^)^ G^''\z) (57) 

for some nontrivial function C^^^ [z) (the nth SRG estimate of the true limiting 
distribution). This result, which can be derived from Eq. (55), explains why 
random fractals, such as critical spanning clusters in percolation, have very 

large mass fluctuations, comparable to the mean, in contrast to normal ran- 
dom objects, such as supercritical spanning clusters, which have "square-root 
fluctuations" (at the scale of the square root of the mean) as a consequence 
of the CLT. 
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Fig. 5. Limiting probability distribution (scaled to have unit mean) for the total 
mass of clusters spanning a square system in one direction (solid line, obtained from 

simulations on a 300 x 300 site lattice with p = 0.592746 ~ Pc) compared with 2x2, 
3x3, and 4x4 SRG approximations (dashed lines, obtained by numerically iterating 
the exact recursions given in the main text). 



In Part III, it is shown that for percolation quantities, the SRG critical fixed 
point, G{z) = limn^Qo G'^"'\z) , is generally not Gaussian. Its shape depends 
sensitively on macroscopic constraints, such as spanning or boundary condi- 
tions, although it is presumably universal for different lattices in the same 
embedding dimension. The fact that shape of the critical order-parameter dis- 
tribution depends on macroscopic constraints, for arbitrarily large systems, is 
a striking signature of long-range correlations, which has apparently not been 
previously noted in the literature. This is also the first conclusion of the FGLT: 
In the central region, the limiting mass distribution of a random unifractal is 
non- Gaussian and depends sensitively on macroscopic constraints. 



Preliminary numerical evidence indicates that satisfactory approximations of 
critical distribution functions can be obtained using remarkably small cells. 

These approximations are produced by numerically iterating the recursions 
given above, either for the probability density directly (which involves discrete 
convolutions) or for one of the generating functions (which involves simple 
multiplications followed by an inverse transform). In Fig. 5, we see that the 
very simple n — 2 x 2 cell scheme discussed above is surprisingly close to the 
results of large-scale simulations for clusters spanning a square site lattice in 
one direction. Even more striking is the good accuracy achieved using only 
slightly larger cells, n = 3 x 3 and n = 4x4. These approximations and others 
are analyzed in more detail in Part III. 



27 



The SRG also predicts universal stretched-exponential decay for the tails of 
the critical limiting distributions, 

-Aiz)z^ as.^oo 
log^(^) ~ <; ' \ (58) 
-B{z)z'^ as z^O 



where g{z) = G'{z) is the limiting probability density function in the limit of 
infinite cell size, and A{z) and B{z) are periodic functions of logz, which are 
essentially constant in practice [13]. (Note that stretched-exponential decay 
of the left tail as 2; ^ holds only if a;~ > 1.) The tail exponents, S and rj, 
can be expressed in terms of other critical exponents and do not depend on 
macroscopic constraints. If the mean scales like /Xn(Pc) ~ and the largest 
possible value scales like ~ n^^ , then the exponent for the right tail is 
given by. 



Similarly, if the smallest possible value scales like x„ ~ n*^ , then the exponent 
for the left tail is given by, 

f, = J—^. (60) 



In the case of spanning clusters in d dimensions, we have S — Df/d, 5'^ — 1, 
and 5~ — 1/d, in which case the tail exponents are. 



and 



(62) 



These conclusions follow directly from classical mathematical theorems on 
Watson-Galton branching processes [33,34,67], as first pointed out by Hovi 
and Aharony [13]. Note that the co-dimension, d — Df, governs the right 
tail, because rare, very large clusters have the dimension of the system d^ 
rather than the normal fractal dimension, Df. Similarly, the co-dimension, 
Df — 1, governs in the left tail, because rare, very small clusters (conditional 
on spanning) are roughly linear objects with dimension one. 
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Subcritical Phase Supercritical Phase 

Fig. 6. Topology of the SRG flow for in the space of scaled probability distributions 
of any structural quantity conditional on spanning, such as the total cluster mass 
or the minimum chemical distance connecting two opposite boundaries. 

These results lead us to the second general conclusion of the FCLT: Outside the 
central region, the limiting distributions of random unifractals have stretched 
exponential tails with universal exponents given by Eqs. (59)-(62). In ongoing 
work [31], the generality of the FCLT is explored by comparisons with perco- 
lation, random graphs, and the Ising model, and the subtle question of what 
happens above the upper critical dimension is also addressed. 



8 Discussion 

Combining the results of this article, we arrive at a fruitful intuitive picture 
of the SRG "flow" (dashed lines) in the space of probability distributions, 
as sketched in Fig. 6. As the system size is increased, the physical manifold 
< p < 1 (sohd hues) is advected from the single-site manifold for = 1 
toward the thermodynamic manifold for A" = oo (interior dotted line). The 
SRG accurately describes the universal flow near the two thermodynamic fixed 
points and the crossover manifold connecting them, far away from the p = 0, 1 
fixed point representing discreteness and non-universality. For a fixed value 
of p 7^ 0,Pc; 1) at small sizes the system trajectory passes near the (unstable) 
critical fixed point where the FCLT holds and then crosses over toward the 
(stable) off-critical fixed point, where the CLT holds, once the system size 
exceeds the correlation length. 

The significance of the SRG is both practical and fundamental. In real ap- 



29 



plications, critical phenomena such as percolation occur in a finite systems, 
which are often small enough compared to the correlation length that finite- 
size effects are important. In such cases, the SRG provides a simple analytical 
or numerical method to approximate the mean, variance, and distribution of 
an order parameter for arbitrary system sizes and values of the control param- 
eter, even when the correlation length and system size are comparable. These 
approximations should be fairly accurate as long as the correlation length and 
the system size are much larger than one. 

Preliminary numerical results for percolation suggest that reasonable accuracy 
can be achieved with remarkably small cell sizes (6 > 2). Developing the best 
small-cell approximations for various quantities in percolation and other sta- 
tistical models will require more extensive comparisons with simulation data 
(as in upcoming parts of this series). Whenever reasonable small-cell schemes 
can be identified, however, they allow trivial calculations which may suffice 
to replace direct simulations in many applications. From a practical point 
of view, simple approximations for a broad range of conditions can be much 
more useful than exact results for certain limiting cases because models like 
percolation are always themselves only crude approximations of real physical 
systems. 

The SRG also provides basic insights into critical phenomena and random 
fractals. The connection with branching processes gives a simple "population- 
growth" picture of how crossover occurs in a continuous phase transition. The 
fixed points of the SRG also explain the nature of fiuctuations in different 
thermodynamic phases, e.g. governed by the CLT and the FCLT. These con- 
clusions have been reached here in the specific context of percolation, which 
is perhaps the simplest case, but it should be possible to extend the SRG to 
other critical phenomena, such as the Ising model, by postulating and ana- 
lyzing appropriate non-stationary branching processes expressing statistical 
self-similarity. 
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